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ABSTRACT 

We report the results of observations of ten rotational transitions of water vapor toward the carbon- 
rich AGB (asymptotic giant branch) star IRC+10216 (CW Leonis), carried out with HerschePs HIFI 
instrument. Each transition was securely detected by means of observations using the dual beam 
switch mode of HIFI. The measured line ratios imply that water vapor is present in the inner outflow 
at small distances (< few x 10 14 cm) from the star, confirming recent results reported by Decin et al. 
from observations with Herschel's PACS and SPIRE instruments. This finding definitively rules out 
the hypothesis that the observed water results from the vaporization of small icy objects in circular 
orbits. The origin of water within the dense C-rich envelope of IRC+10216 remains poorly understood. 
We derive upper limits on the H^O/H^O and H^O/H^O isotopic abundance ratios of ~ 5 x 10~ 3 
(3 <r), providing additional constraints on models for the origin of the water vapor in IRC+10216. 
Subject headings: circumstellar matter - stars: AGB and post- AGB - stars: abundances 



1. INTRODUCTION 

The carbon-rich AGB (asymptotic giant branch) star 
IRC+10216 (C W Leonis) is surrounded by a dense, out- 
flowing envelope that harbors a rich chemistry. Be- 
cause of its proximity to Earth (D ~ 100 — 170 pc) and 
its large mass-loss rate (~ 1.5 — 3 x 10~ 5 M Q yr _1 ; see 
Crosas & Menten 1997 and Schoier & Olofsson 2000), 
IRC+10216 has proven to be the most valuable source 
available for the study of chemistry in a carbon-rich as- 
trophysical environment; more than 60 molecules (e.g. 
He et al. 2008) have been detected within its circum- 
stellar envelope (CSE), made up of the elements H, C, 
N, O, F, Na, Mg, Al, Si, S, P, F, and K. The photo- 
sphere of IRC+10216 has been significantly enriched in 
carbon - owing to the dredge-up of helium-burning nu- 
clear reaction products from the core - leading to a large 
C/O ratio (~ 1.4; Winters, Dominik, & Sedlmayr 1994) 
and giving rise to a circumstellar chemistry that is qual- 
itatively different from that typical of oxygen-rich AGB 
stars. In particular, while the most abundant molecules 
in O-rich CSEs are CO and H 2 0, those in IRC+10216 
are CO, HCN, and C 2 H 2 . Prior to 2001, the only oxygen- 
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rich molecules reported in the CSE of IRC+10216 were 
CO, SiO, and HCO+ (Lucas & Guelin 1999), but in the 
last decade four other O-bearing molecules have been 
detected: H 2 (Melnick et al. 2001), OH (Ford et al. 
2003), H 2 CO (Ford et al. 2004), and C 3 (Tenenbaum 
et al. 2006). 

The discovery of water vapor in IRC+10216, obtained 
by observations of the lio — lot rotational transition near 
557 GHz with the use of the Submillimeter Wave As- 
tronomy Satellite (SWAS), was puzzling because mod- 
els for thermochemical equilibrium in the photospheres 
of carbon-rich stars predict water abundances of only 
~ 10~ 10 relative to H 2 (e.g. Cherchneff 2006), many or- 
ders of magnitude below the detection threshold or the 
abundance ~ 10 -7 inferred (Agundez and Cernicharo 
2006; Gonzalez- Alfonso, Neufeld and Melnick 2007) from 
SWAS observations. Several hypotheses have been pro- 
posed for the origin of the observed water vapor. Ford 
& Neufeld (2001) proposed a model whereby the water 
vapor was released into the outflow by the vaporization 
of a Kuiper belt analog, in which orbiting icy objects - 
heated by the enhanced luminosity of the AGB star - un- 
derwent sublimation. Willacy (2004) subsequently sug- 
gested that the observed water could be formed from CO 
and H 2 by means of Fischer- Tropsch catalysis on metal- 
lic dust grains, whilst Agundez and Cernicharo (2006; 
hereafter AC06) proposed that water might be formed 
in the outer envelope by means of radiative association 
of H 2 with O atoms produced by the photodissociation 
of 13 CO exposed to the ultraviolet interstellar radiation 
field (ISRF). Because SWAS had access to only a single 
transition of water vapor, it was unable to place observa- 
tional constraints that could readily distinguish between 
these hypotheses. 

In a theoretical study (Gonzalez- Alfonso, Neufeld and 
Melnick 2007; hereafter GNM) performed in anticipation 
of the Herschel Space Observatory, it was discussed how 
each of these hypotheses gave rise to a specific prediction 
for the spatial distribution of the emitting water vapor. 
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GNM investigated how multitransition observations with 
Herschel might be used to determine that distribution, 
and showed, in particular, that the relative strength of 
lines of higher excitation than those accessible to SWAS 
was a decreasing function of the inner radius, R ln , of the 
region containing water vapor. The vaporization of a 
Kuiper belt analog would lead to i?; n ~ 2 x 10 15 cm ~ 
30i?» (Model B in GNM; where i?» - 8 x 10 13 cm is the 
assumed stellar radius) because any icy object (at least 
of small size and on a circular orbit) within that radius 
would have been vaporized before IRC+10216 ascended 
the AGB. Water production by means of Fischer- Tropsch 
catalysis would result in a similar value of R m , accord- 
ing to calculations of Willacy (2004). The formation of 
water following the photodissociation of CO in an outer 
shell, as suggested by AC06, by contrast, would lead to 
a much larger R in ~ 4.3 x 10 16 cm - 500 i?* (Model 
C in GNM), and considerably smaller strengths for the 
higher-excitation water transitions. Conversely, were Ri n 
significantly smaller than 2 x 10 15 cm - a possibility that 
had not been anticipated by any specific model - then the 
high-excitation transitions would be relatively stronger 
(Model A in GNM). 

In this Letter, we report the results obtained from Her- 
schel (Pilbratt et al. 2010) observations of 10 rotational 
transitions of water vapor, carried with the Heterodyne 
Instrument for the Far Infrared (HIFI; de Graauw et 
al. 2010). The observations and data reduction are de- 
scribed in §2, and the spectral line profiles and line in- 
tensities are presented in §3. In §4, we discuss the spa- 
tial distribution inferred for the water vapor in the CSE 
of IRC+10216, in the context of various hypotheses for 
its origin. We compare our results with those reported 
recently in an entirely independent study performed by 
Decin et al. (2010; hereafter D10) with the use of the 
PACS and SPIRE instruments on Herschel. 



2. OBSERVATIONS AND DATA REDUCTION 

Observations of IRC+10216 were carried out in May 
2010 as part of the HIFISTARS Key Program. We 
used the HIFI instrument in dual beam switch (DBS) 
mode to target 10 rotational transitions of water va- 
por. The list of observed transitions, and the details 
of each observation, are given in Table 1. The tele- 
scope beam was centered on IRC+10216 at coordinates 
a = 9h 47m 57.38s, 5 = +13°16'43.7" (J2000), and the 
reference positions were located at offsets of 3' on ci- 
ther side of the source. The data were processed using 
the standard HIFI pipeline to Level 2, providing fully 
calibrated spectra with the intensities expressed as an- 
tenna temperature and the frequencies in the frame of 
the Local Standard of Rest (LSR). For five of the ten 
observed transitions, we used a pair of local oscillator 
(LO) frequencies, separated by a small offset, to confirm 
the assignment of any observed spectral feature to either 
the upper or lower sideband of the (double side band) 
HIFI receivers. The resultant spectra were coadded so 
as to recover the signal-to-noise ratio that would have 
been obtained at a single LO setting. Spectra obtained 
for the horizontal and vertical polarizations were found 
to be very similar in their appearance and noise charac- 
teristics and were likewise coadded. 
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Fig. 1. — Hcrschcl/HIFI spectra of six rotational transitions of 
H2O with upper state energies Eu/k < 160 K. The blue line super- 
posed on each spectrum shows the profile of the lio — loi ground- 
state transition of ortho- water (top left panel). Doppler velocities 
are expressed relative to the LSR velocity of the source, taken as 
-25.5 kms" 1 . 
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Fig. 2. — same as Fig. 1, except for four rotational transitions of 
H2O with upper state energies E\j/k > 160 K, and one rotational 
transition each of CO, 13 CO, SiO and SiS. 
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3. RESULTS 

The integrated antenna temperatures and line fluxes 
are given for each detected transition of water vapor in 
Table 2, together with the energy of the upper state. The 
HIFI instrument used for these observations provides a 
spectral resolution (1.1 MHz for the Wide Band Spec- 
trometer) that is much superior to that of the PACS or 
SPIRE instruments; this property of HIFI typically pro- 
vides for the unambiguous identification of strong lines, 
and yields spectral line profiles that can be used to probe 
the kinematics of the outflowing gas. In Figures 1 and 2, 
we present the spectra obtained for each of the ten de- 
tected water vapor transitions, each smoothed to a spec- 
tral resolution of 1 kms -1 to improve the signal-to-noise 
ratio. The lower four panels in Figure 2 show compari- 
son spectra for one transition each of the molecules CO, 
13 CO, SiO and SiS; taking advantage of the wide band- 
width of the HIFI receivers (4 GHz for v < 1250 GHz 
and 2.4 GHz for v > 1410 GHz), we observed these lines 
simultaneously with nearby water transitions. In Fig- 
ures 1 and 2, the blue line superposed on each spectrum 
shows the profile of the lio — loi ground-state transition 
of ortho-water. 

Although a full analysis of the spectral line profiles will 
be the subject of a future publication, several features are 
obvious from an inspection of the spectra. Most of the 
water lines share a distinctive profile that is quite dissim- 
ilar from those of either the 12 CO or 13 CO transitions. In 
particular, the blue sides of the water lines are very sim- 
ilar to the nearly parabolic profile shown by 12 CO, while 
the red sides are more similar to that of the 13 CO line, 
which shows the more rectangular profile expected of an 
optically-thin line. In these respects, the typical water 
line profiles are most similar to that of the silicon monox- 
ide line shown in Figure 2. The observed asymmetry in 
the line profiles likely reflects an asymmetry in the distri- 
bution of the emitting H2O and SiO, perhaps suggesting 
a common origirP^I. The highest excitation water line to 
have been detected in our study, the 1162.912 GHz tran- 
sition, exhibits a profile that is rather different from the 
typical water line profile, with a narrow core that is sim- 
ilar to that present in the SiS 1102.029 GHz line profile; 
here, emission from the acceleration zone may be impli- 
cated. In addition, the other water lines typically show a 
small narrow emission bump near the systemic velocity 
of the source, perhaps also representing emission from 
material that has not yet been fully accelerated. 

4. DISCUSSION 

In Figure 3, we compare the measured line fluxes tabu- 
lated in Table 2 with the predictions of the GNM models. 
Here, the observed line fluxes are represented by black 
crosses for each of the ten detected transitions, ordered 
from left to right in increasing energy of the upper state, 
and the GNM predictions are shown by the solid lines 
for various values of the inner radius, i?i n , of the water 
emitting region. These predictions were "calibrated" by 

11 Cherchncff (2006) has proposed pulsation-driven shock waves 
as the source of SiO in carbon-rich AGB stars, and - although the 
models presented in that 2006 study failed to predict any enduring 
enhancement in the water abundance behind such shock waves - 
preliminary results from a more recent calculation do indeed sug- 
gest (Cherchncff 2010) that the observed water vapor could be 
produced along with SiO by shocks in the inner envelope. 



TABLE 1 
Observations of IRC+10216 



Transition 


V 

(GHz) 


Mixer 
band 


HPBW a 

(") 


Date 
2010 


t b 

obs 

(s) 


lio 


- 1()1 (o) c 


556.936 


lb 


38 


May 4, 11 


3150 d 


2n 


- 2 02 (p) 


752.033 


2b 


28 


May 12 


3236 d 


2fj2 


- Ill (p) 


987.927 


la 


22 


May 16 


3242 d 


3l2 


- 3 3 (0) 


1097.365 


4b 


19 


May 11 


10446 d 


111 


- Ooo (p) 


1113.343 


4b 


19 


May 11 


3200 d 


3l2 


- 2 2 i (0) 


1153.127 


5a 


18 


May 12 


1539 e 


3 2i 


- 3l2 (0) 


1162.912 


5a 


18 


May 12 


1539 e 


221 


- 2i2 (0) 


1661.008 


6b 


13 


May 14 


1615-f 


2l2 


- loi (0) 


1669.905 


6b 


13 


May 14 


1615-f 


303 


- 2l2 (0) 


1716.770 


7a 


12 


May 14,15 


10412 d 



a Half power beam width 

b Total observing time, including overheads 

c The letters o and p indicate whether the transition is of ortho- 
or of para-water 

d Data obtained with two separate LO settings of duration t b s /2 
e,f Pairs of lines observed simultaneously 

adjusting the assumed water abundance to match the 557 
GHz line flux detected by SWAS. Clearly, the best fit to 
the data is obtained for the smallest value of R m con- 
sidered by GNM, 4.5 x 10 14 cm (~ 5.6 R„), but even for 
that value the model underpredicts the relative strengths 
of the transitions of highest excitation. Thus, our HIFI 
observations confirm the inference drawn by D10 - from 
the detection of high-lying water rotational lines with 
PACS and SPIRE - about the spatial distribution of wa- 
ter vapor; water is clearly present at distances smaller 
than 4.5 x 10 14 cm (~ 5.6 i?*) from the star. The pres- 
ence of water that close to the star definitively falsifies 
the model proposed by Ford & Neufeld (2001), in which 
the origin of the observed water vapor was the vaporiza- 
tion of icy comets in a Kuiper belt analog, because any 
icy object (at least of small size and on a circular orbit) 
within 4.5 x 10 14 cm would have been vaporized before 
IRC+10216 ascended the AGB. The models proposed by 
Willacy (2004; i.e. Fischer- Tropsch catalysis) and par- 
ticularly AC06 (production via radiative association in 
an outer shell) are similarly excluded. The fact that the 
1113 GHz / 557 GHz line ratio is in good agreement 
with the GNM model implies that the ortho-to-para ra- 
tio is close to 3, the value assumed in the model. The 
dashed black line in Figure 3 shows our best-fit model, 
with i?; n = 1.0 x 10 14 cm (~ 1.3 i?*), and a slightly larger 
water abundance H2O/H2 = 8.1 x 10 -8 than that as- 
sumed in GNM, now chosen to optimize the fit to all the 
transitions observed with HIFI. 

At present, the origin of water vapor in IRC+10216 
remains poorly-understood. D10 have proposed an alter- 
native model - discussed in greater detail by Agundcz, 
Cernicharo and Guelin (2010; hereafter ACG) - involving 
photochemistry in the inner envelope. As in the AC06 
model, oxygen atoms are generated by photodissociation 
of 13 CO and SiO by the ultraviolet ISRF, but - unlike 
the AC06 model - the UV radiation is assumed to pene- 
trate deeply into a clumpy CSE. In the D10/ACG model, 
the oxygen atoms are liberated close to the star, where 
the temperature is sufficient to drive H2O production via 
a sequence of two H-atom extraction reactions with ac- 
tivation energy barriers: 0(H2, H)OH(H2, H)H20. Be- 
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Fig. 3. — Comparison of the measured water line fluxes (black 
crosses) with predictions of the GNM models for various values 
of the inner radius, R ln , of the water-emitting region. The ten 
rotational transitions appear from left to right in order of increasing 
energy for the upper state (except for the 1153 GHz and 1097 GHz 
transitions which originate in the same upper state). 



TABLE 2 

Water line fluxes measured toward IRC+10216 



Transition 


V 


E v /k 


fT A dv 


Flux a ' 6 




(GHz) 


(K) 


(Kkms- 1 ) 


(10-2° W cm" 2 ) 


lio - loi (o) c 


556.936 


26.7 


10.3 


0.89 


2n - 2 2 (p) 


752.033 


136.9 


6.0 


0.70 


2()2 - 111 (p) 


987.927 


100.8 


7.2 


1.11 


3l2 - 303 (o) 


1097.365 


215.2 


10.7 


1.83 


In - Ooo (p) 


1113.343 


53.4 


14.6 


2.53 


3l2 - 221 (o) 


1153.127 


215.2 


3.4 


0.62 


321 - 3l2 (o) 


1162.912 


271.0 


5.3 


0.96 


221 - 2i2 (o) 


1661.008 


159.9 


15.5 


4.09 


2l2 - loi (o) 


1669.905 


80.1 


35.1 


9.56 


303 - 2l2 (o) 


1716.770 


162.5 


24.2 


6.74 



a For an unresolved source at the beam center 

b We conservatively estimate the flux uncertainty to be less than 
15%. 

c The letters o and p indicate whether the transition is of ortho- 
or of para-water 



cause the mean visual extinction through the CSE is 
^100 mag, this scenario requires the existence of chan- 
nels of greatly reduced extinction through which the ul- 
traviolet radiation can penetrate. The consequences of 
this model for the abundances and spatial distribution of 
the many other molecules detected in IRC+10216, many 
of which have been observed interferometrically, has yet 
to be fully investigated. 

One possible test of this model might involve a search 
for H^O or H^O. Because the 12 CO photodissociation 
rate is sharply reduced by self-shielding, D10 have em- 
phasized the importance of 13 CO as a source of atomic O 
that can react to form H2O. The photodissociation rates 
for C 17 and C 18 would presumably be at least as large 
as that for 13 CO. Thus, if 13 CO, C 17 and C 18 were 
the only suppliers of atomic oxygen, the H^O/H^O 
and H^O/H^O ratios would approach the 13 CO/C 18 
and 13 CO/C 17 ratios respectively. Given the 13 C/ 12 C, 
18 0/ 16 and 17 0/ 16 isotopic ratios determined by Ka- 
hane et al. (1992) and by Cernicharo, Guelin & Kahanc 
(2000) for IRC+10216, the 13 CO/C 17 and 13 CO/C 18 
ratios are respectively ~ 18 and 28, each a factor of 45 
(= 12 C/ 13 C) smaller than the elemental 16 0/ 17 and 
16q/18q ratios in the CSK jf SiQ or 12 C0 arc significant 

additional sources of atomic oxygen in the D10/ACG 
picture, then a H^O/H^O ratio larger than 18 - or a 
H^O/H^O ratio larger than 28 - might still be consis- 
tent with the model. 

Data obtained in a full HIFI spectral survey carried out 
toward IRC+10216 (Cernicharo et al. 2010) place upper 
limits on the flux of the H^O (552.021 GHz) and H* 8 
(547.676 GHz) lio — loi transitions. Comparing these 
with the observed flux in the 556.936 GHz H* 6 li - l i 
transition, we determined that the spectral survey places 
3 a lower limits on both the 556.936 GHz/552.021 GHz 
and 556.936 GHz/547.676 GHz line ratios of - 100. 
Taking account of optical depth effects, we find that 
these correspond to a lower limit of ~ 200 on both the 
H^ 6 0/Hf O and H^O/H^O abundance ratios. Given 
the elemental isotopic ratios 16 0/ 18 ~ 1260 and 
16 0/ 17 - 840 (Kahane et al. 1992), these limits imply 
that the abundance of the minor isotopologucs could only 
be enhanced (by means of isotope-selective photodissoci- 
ation of CO, for example) by at most a factor 4 (H^O) 
or 6 (H3, 8 0). Detailed modeling will be required to de- 
termine whether our non-detections of H^O and H^O 
are consistent with the water production mechanism pro- 
posed by D10/ACG. 

HIFI has been designed and built by a consortium of in- 
stitutes and university departments from across Europe, 
Canada and the United States under the leadership of 
SRON Netherlands Institute for Space Research, Gronin- 
gen, The Netherlands and with major contributions 
from Germany, France and the US. Consortium mem- 
bers are: Canada: CSA, U. Waterloo; France: CESR, 
LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, 
MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, 
Osservatorio Astrofisico di Arcctri- INAF; Netherlands: 
SRON, TUD; Poland: CAMK, CBK; Spain: Observato- 
rio Astronomico Nacional (IGN), Centro de Astrobiologa 
(CSIC-INTA). Sweden: Chalmers University of Tech- 
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nology - MC2, RSS & GARD; Onsala Space Observa- 
tory; Swedish National Space Board, Stockholm Univer- 
sity - Stockholm Observatory; Switzerland: ETH Zurich, 
FHNW; USA: Caltech, JPL, NHSC. 
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